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Abstract 

We explore various pitfalls and challenges in determining the equation- 
of-state (w) of dark energy component that dominates the universe 
and causes the current accelerated expansion. We demonstrated in an 
earlier paper the existence of a degeneracy that makes it impossible to 
resolve well the value of w or its time-derivative with supernovae data. 

Here we consider standard practices, such as assuming priors that w 
is constant or greater than —1, and show that they also can lead to 
gross errors in estimating the true equation-of-state. We further con¬ 
sider combining measurements of the cosmic microwave background 
anisotropy and the Alcock-Paczynski test with supernovae data and 
find that the improvement in resolving the time-derivative of w is 
marginal, although the combination can constrain its present value 
perhaps to 20 percent uncertainty. 

PACS number(s): 98.62.Py, 98.80.Es, 98.80.-k 
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1 Introduction 


Measurements of Type lA supernovae have shown that the expansion of the 
universe is accelerating ||l|, |^, suggesting that most of the energy density of 
the universe consists of some form of dark energy with negative pressure [^]. 
Combining measurements of the cosmic microwave background anisotropy 
and observations of large-scale structure provides important corroborating 
evidence. ii Two candidates for the dark energy are a cosmological con¬ 
stant (or vacuum density) and quintessence,]^] a time-varying, spatially in¬ 
homogeneous component. In a previous paper [0 (Paper I), we addressed 
the question of whether supernova measurements can be used to measure the 
equation of state (EOS) of the negative pressure component, the ratio w of 
the pressure to the energy density. The issue is important because w = —1 
for a cosmological constant whereas w takes on different values and can be 
significantly time-varying in the case of quintessence. P, |] Under the assump¬ 
tion that w is constant, its value can be determined to better than 5 per cent 
by measuring several thousand supernovae distributed equally between red 
shift z = 0 and z = 2. However, we showed that a degeneracy opens up if 
w is time-dependent which makes it impossible to determine accurately the 
current value of w or its time-derivative. The cause of the degeneracy is that 
supernovae measure luminosity distance, which is related by a multi-integral 
expression to the EOS as a function of red shift, w{z). Widely different w{z) 
can have the same multi-integral value. 

The purpose of this paper is to explore some pitfalls and challenges in 
determining the EOS and its time variation using supernovae. For example, 
we shall show how the standard practice of considering only models with 
w > —1 or only models with constant w when doing likelihood analyses can 
lead to grossly incorrect results. For example, we will illustrate cases where 
the standard practice will suggest that w is near -1 or much more negative 
than —1 when, in fact, w is significantly greater than -1 and rapidly time 
varying. We shall show that a non-zero value of dw/dz is more easily de¬ 
tected if dw/dz > 0 than if dw/dz < 0. We shall also contrast measuring w 
for the negative pressure component alone {wq) versus the mean value for the 
total energy density (including ordinary and dark matter), wt- The degen¬ 
eracy problem is less severe for wt, but this parameter provides less useful 
information. We consider possibilities of breaking the degeneracy between 
w and dw/dz by combining supernovae results with either cosmic microwave 
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background anisotropy measurements and/or the Alcock-Paczynski test. We 
shall show that neither additional test signihcantly improves the measure¬ 
ment of the time-variation of w, although optimistic assumptions about the 
Alcock-Paczynski test suggest that the current value of w can be measured 
to within 20 percent or so. 

We conclude that a new, yet to be found test has to be devised to resolve 
well the cosmic EOS and its time variation. We stress that with current 
data it is possible to determine the EOS to about a factor of two. For a 
future experiment to signihcantly enhance the determination of the EOS, 
and enable the distinction between a constant EOS and a time-dependent 
one, it needs to resolve the equation of state at the 10% level or better. 

The results of our analysis agree with many other analyses |^, 


T^, |T^, |I^ [T^, [T^, ^ , although not always with their interpretation. 


and can be used to explain why some other analyses seem to indicate a 
superior resolving power of SN measurements alone ^ 
combination with other measurements 
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Some of the latter 

analyses implicitly assume unrealistic accuracy in independent determination 
of cosmological parameters, by not including self-consistently the uncertainty 
in w{z) in all measurements. For example, assume a reported resolution of 
matter energy density which was based on assuming wq = — 1. 

Type lA supernovae have intrinsic variability of about 0.15 in absolute 
magnitude, but currently the errors in measuring distant SN are above this. 
There are ongoing programs to extend the search to deeper redshifts and 
improve measurement quality. The proposed SNAP satellite plans to measure 
2000 SN per year, mostly in the range of redshifts 0.1 < < 1.2, and 
some as far as z = 1.7. The anticipated error in individual magnitudes 

is Am = 0.15 statistical and 0.02 systematic error, which yields about 1% 
relative error in luminosity distance di. 

For our numerical estimates, we have generated 50 SN magnitudes ran¬ 
domly chosen from a uniform distribution in z values, between z = 0.1 and 
z = 2. Magnitudes were generated from a Gaussian distribution with mean 
value m{z) calculated using hducial models. We have used our 50 points to 
simulate approximately 2000 SN by reducing their magnitude error by a fac¬ 
tor a/40 to 0.03 from the minimal 0.15 magnitude. Thus each generated point 
corresponds to 40 SNAP-like points, binned together. This corresponds to 
1.4% relative error in Hence, our analysis is based on a SN search more 
extensive than the actual SNAP proposal. To obtain a quantitative estimate 
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of how well models are resolved, we use one of two procedures. First, we can 
find the maximum likelihood contours of the various models for each of the 
hducial models and explore the degeneracy in parameter space. Alternatively, 
we can assume that all models which predict dL^z) within 1% of the hducial 
cosmological model for all between 0 and 2 are deemed indistinguishable. 
We hnd that both approaches give comparable results. That is, the 95% CL 
likelihood contours using the hrst procedure are roughly equivalent to the 
indistinguishability region of the second. 


2 Dependence of luminosity distance on dark 
and total EOS 


Luminosity distance is dehned to be the ratio of luminosity C to hux JF 


di = 



(1 + z)r, 


( 1 ) 


where the present value of the scale factor oq is normalized to unity (through¬ 
out subscript 0 denotes present values); r is the coordinate distance 

dx 



and H is the Hubble parameter H = a/a. The observed SN magnitudes are 
related to di, 

m{z) = M + 25 + 5 log^o [HodL{z )], (3) 

M being the SN absolute magnitude. 

For a hat universe with two energy sources, matter (including dark mat¬ 
ter) and a dark Q-component, there are two equivalent routes to computing 
di without assuming any prior about the time dependence of wq. One way is 
to use the algebraic relation between the total energy density and the Hubble 
parameter H. Using conservation equations, the energy densities of the dark 
component pq and that of ordinary matter are given by 


= (P( 3 )o(l + zf exp 
= (Pm)o(l + z)^. 
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where 5 ^ denotes the present ratio of matter to dark energy densities g = 
(|^)o, H can be expressed as 

1/2 

(7) 


H = Ho{l + 
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Substituting this into gives 

l + z 


dx 

3 / wq — 


X 


dr, = 


(1 + ^) 


[1 + 9?'^ j 


dx 

3/2 


a; 


g + exp 3 


1 / 


1 -1/2 


( 8 ) 


An equivalent approach is to treat the sum of dark matter and the dark 
energy component as a single cosmic fluid with average equation of state 
wt{z) = Q.q{z)wq{z), where 

Wq 


Wt{z) = 


1 + {Pm/P q) 


Wq 


l + gexp -3//+^'«;q(x) 


> dx 


( 9 ) 


Since is proportional to the total energy density in the universe, it can 
be expressed in terms of wt as follows 
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— = exp 

rlnj 


/• 1 +^ _ . ..dx 
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Here we have used the conservation equation for the total energy density. 
Using ([10|) we can express dz, in terms of wq as follows, 
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This expression for rf/, as a function of wq has one more integral than the 
relation in Eq. (|[), but it is pedagogically useful in demonstrating that is 
sensitive only to a weighted average of wt or wq and not to their detailed 
time dependence. 


3 Constraining dark and total EOS using SN 
measurements 

Based on the previous sections, a number of lessons can be learned about 
measuring the EOS. First, the relation between wt and in ([TT| ) involves 
an integral, so we do expect some degeneracy in the determination of wt{z) 
from SN measurements. To determine the EOS of the dark energy itself, wq, 
the total energy density must be resolved into a matter component and a 
dark energy component. Hence, if flm and Oq are not known from indepen¬ 
dent measurements, determining wq entails an additional uncertainty. For 
example, consider a flat universe with Oj’ = = 1- Since the matter 

EOS is Wm = 0, it follows that wt = wqQq. Two models with different 
values of wq may produce the same value of wt and, consequently, d^, due 
to offsetting differences in the value of flm- 

The more negative wt is, the faster is the expansion. Therefore, a more 
negative (positive) wt will make di larger (smaller). In addition, light from 
earlier times (emitted at higher z values) must pass through the low z uni¬ 
verse to reach us. This means that changes in wt at lower z affect d^ at 
higher z. The converse is not true. Changes in wt at high z do not affect d^ 
at lower 2 ;. 

Consequently, it is not surprising that SN measurements of d^ provide 
stronger constraints on wt{z) at low z than at high 2 ;. In particular, if all 
cosmic parameters other than wq{z) are hxed, there is a particular, relatively 
low value of z = z* for which wt{z) is most tightly constrained. This value 
of z* is clearly seen in our numerical results and was noted independently 
by [^, |26|. For example. Figure |l| shows the EOS wt{z) for three models 
each of which is obtained by best £t to d^i^z) for a fiducial model {wq, Qm) = 
(—0.7 — O.Sz, .3) using one of three fitting assumptions: (1) that the dark 
EOS is constant; (2) that the dark EOS is linear wq = Wq+Wiz; and, (3) that 
the total EOS is quadratic wt = A + Bz -f Cz^. While the real degeneracy 
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red shift z 


Figure 1: wt{z) for three best fit models of three hts under three different 
assumptions: constant wq (dashed), linear wq (solid), and quadratic wt 
(dotted), to data generated from a single hducial model; {wq, = (—0.7 — 
O.Sz, .3). All hts prefer = wt(z* ~ .15) ~ —.52, but diverge for other 
values of 2 ;. 
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is stronger than what is seen in the hgure, we have chosen three examples 
to illustrate the existence of and z*. As can be seen, the hts disagree 
signihcantly for z far from z* = 0.15, but all hts agree near z*. 

Unfortunately, the resolution of wq{z*), the quantity which most inter¬ 
ests us, is degraded when we do not £x flrn but, instead, allow for the current 
uncertainty in its value. In Fig. Q, we show some linear hts to simulated data 
generated from the hducial model {wQ,Qrn) = (—I,-3). The hts are repre¬ 
sentative examples which ht the hducial model to within the 95% confidence 
region. The upper plot shows that wt{z) (with flm hxed at 0.3) is relatively 
well resolved, and particularly well resolved at around redshift 2 :* = .3. The 
resolution is not that sharp in the middle plot which shows the corresponding 
wq, but a special point of enhanced resolution around z = A is still clearly 
seen. If one lets vary in the realistic range of 0.2 — 0.4, then wq becomes 
poorly resolved and the spread at z* increases signihcantly. Similarly the 
spread in w* and z* increases signihcantly if more general functional forms 
of the EOS are considered. 

We would like to stress that the constant or linear forms of wq that we 
use are not meant to be anything more than simple concrete examples to 
highlight the fact that we are dealing with a degenerate parameter space. 
Showing that if one assumes a linear wq{z), then it can be resolved to, 
say, 50% does not logically mean that it can be measured to 50% accuracy 
generally since wq{z) is resolved with diherent accuracy depending on its 
functional form. This can be illustrated with the following examples. In 
Fig. 1^, the diherence in magnitude (Am) for models with various EOS is 
shown. There are three clusters of points, each of which corresponds to a 
simulation of SN data for pair of diherent models. Each pair consists of a 
constant and linear wq. Each pair can be clearly separated from other pairs 
but the constant and linear “members” of a pair cannot be distinguished by 
SN data. The examples chosen for Fig. ^ have unrealistic large derivatives 
(of order unity) and therefore start to diverge from their constant partners 
for large z. More realistic examples with smaller derivatives or oscillatory 
behaviour will be much harder to distinguish from a constant EOS. 

Clearly the treatment of the SN analysis is important. If it is assumed 
that Wq is constant, the hgure shows that diherent values can be resolved 
to high accuracy, but if the assumption of constancy is relaxed and a linear 
2 : dependence is allowed it becomes clear that the data can determine well 
only a single relation between wq and wi and that wq{z) is poorly resolved. 



.5 1.5 

red shift z 


Figure 2: Models within 95% CL region of a fit to data generated from the 
hducial model {wq, flm) = (—1, 0.3) assuming wq = wq + wiz. Top: The 
total EOS wt{z), for three different linear models. Middle: wq(z) assuming 
Qm = 0.3 exactly, for the same linear models. Bottom: wq{z) for nine 
models, assuming that 0.2 < < 0.4 (no relation between the dashed, 

dotted and solid lines of the bottom g^anel to those of the middle and top 
ones). 






red shift z 


Figure 3: Magnitude differences between pairs of degenerate models and a 
flat pure matter = 1) Universe. Each pair consists of simulated data 
points generated from one constant wq model (open circles) and one linear 
wq model with a large (positive or negative) derivative (full squares). The 
pairs are well separated but it is hard to separate between “members” of 
each pair. 
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Figure 4: 95% CL contours of fits to data generated from two fiducial 

models. The curvature of degeneracy contours is positive for a positive wq 
fiducial model (right) and negative for a negative wq hducial model (left). 

The degree of degeneracy exhibited in the wq = const, hts depends on 
whether wq is positive or negative. Recall that if different models yield a 
total EOS wt = wqQq that is approximately equal, they are degenerate, and 
therefore changes in wq can be compensated by changes in Qq (or equiva¬ 
lently, in Om). The difference between the case where wq is positive is due 
to the specihc way in which this compensation mechanism operates. If wq is 
positive, the curvature of degeneracy lines in {wQ,Qrn) plane is positive, as 
shown in the right panel of Fig. ^ for a fiducial model with wq = -1-0.5. Con¬ 
versely, if Wq is negative, the curvature of the degeneracy line in {wq, Qm) 
plane is negative, as demonstrated in the left panel of figure We have 
found that this result is unaffected by the value of the derivative of wq, even 
if it is quiet large. 
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4 Common Practices and Pitfalls in Deter¬ 
mining wq(z) 

The previous section (and Paper I) show that the determination of wq{z) 
from SN data is a more delicate process than it would seem. If we know a 
priori that wq is constant, then its value can be determined quite accurately. 
However, without this assumption, wq is poorly determined, and matters 
much get much worse if Vtm is uncertain. 

The analysis can be further confounded if certain common practices are 
followed. For example, many analyses assume that wq is constant and pre¬ 
sume that, even if wq is time-varying, the constant- wq £t will provide the 
mean value over recent epochs. Another common practice is to impose the 
condition that wq{z) be limited to — 1 < wq{z) < 1, based on the positivity 
and stability conditions that apply to most (but not all) forms of dark en¬ 
ergy. We shall see that both practices can produce enormous distortions of 
the likelihood surface that lead to grossly incorrect conclusions. 

For example, we have tried to fit data generated from a hducial model 
with Wq = —.7 + .8z, and flm = -3 over a redshift range 0 < z < 2. Note that 
the hducial model has wq > —1 for all z. Yet, if we do a best-£t assuming the 
prior that wq is constant, we find it to be {wq, Qm) — (—1-75, .65). Not only 
does the best-£t have wq < — 1, but the whole 95% conhdence contour lies in 
a region where wq < — 1. The results are the elongated contours in the lower 
part of Fig. The reason for such a strange result can be understood from 
the functional dependence of and wq. Assuming that wq(z) = Wq + Wiz, 
the energy density of the dark component is given in eq.(^, 

PQ = (pQ)o(l + zf^°-^^+^^expl3wiz], (12) 

so increasing wi and decreasing tco have opposite and compensating effects, 
which tend to cancel each other’s inhuence on The exponential factor 
determines the quantitative details of this compensation since changes in wq 
need to compensate also for changes in the exponential factor. It is therefore 
clear that the compensation cannot be perfect over a range of 2;’s. If we insist 
on the prior of constant EOS (he., wi = 0), the htting procedure will pick 
out a value of wq which is much more negative than the hducial value. In the 
figure we have picked a hducial with a large positive derivative to illustrate 
our point, but it is clear from our discussion that the same problem arises 
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when time dependence is weaker, or in cases that the EOS has a more general 
functional form. 

Introducing a prior that wq > — 1 can give a very misleading impression 
of how well Wq is resolved. For example, suppose that we assume the priors 
that Wq is constant and wq > —1, as is standard practice. The results are 
shown in Fig. the small contours truncated at wq = —1. They seem to 
suggest that the data supports the conclusion that wq = — 1 with a high level 
of conhdence. The best £t is {wQ,flrn) — (—1,.58). Yet, this is not related 
in any obvious way to the fiducial model, wq = —0.7 + 0.8z and flm = 0.3. 
The values of per degree of freedom for the best £t models of Fig. are 
reasonable, .95 for the unconstrained £t, and 1.39 for the constrained fit. 
So, what appears to be a compelling result is actually a total distortion. Of 
course, it is also conceivable that the actual wq{z) is less than —1, in which 
case the same procedure of introducing a prior would falsely suggest that 
w = —1 fits well. 


5 Asymmetry in determination of the time 
dependence of the dark EOS 

An EOS in which w{z) has a large positive time derivative is easier to detect 
than one which has a large negative time derivative. In either case, the 
derivative must be large to be detected, as pointed out in Paper I, but here we 
are demonstrating that the challenge is asymmetric. The point is illustrated 
in Fig. ^ The middle panel shows a fiducial model with a modest value 
of wi = 0.2, and, as can be seen, this case cannot be distinguished from 
a model in which wi = 0. That is, the 95% CL contours overlap the line 
tci = 0, corresponding to no time variation. The left and right panels show 
cases in which wi = 0.8 and —0.8, respectively. The contours for the wi = 0.8 
case (left) lie far from the wi = 0 line, so the time variation is detectable 
in this example. On the other hand, the contours for the wi = —0.8 case 
(right) overlay the wi = 0, so the time-variation is not resolved. 

This effect can be explained by considering the variation of the total 
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Figure 5: Constrained (small) and unconstrained (larger and more negative) 
68% and 95% confidence contours of a fit to data generated from a fiducial 
model with linear wq (wQ,Qm) = (—.7 + .8z, .3). The fit is done under the 
(wrong) assumption that Wq is constant. 
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Figure 6: Likelihood contours (68% (lighter) and 95% (darker) C.L.) in 
the {wo, wi) plane, for hts to data generated from 3 different hducial mod¬ 
els. LEFT: = (—0.7,0.8,0.3). MIDDLE: {wo,Wi,Qrn) = 

(—0.7, 0.2, 0.3). RIGHT: (wo,Wi,ilrn) = (—0.7,—0.8, 0.3). Only the results 
shown in the left panel are inconsistent with a constant (rci = 0) wq model. 


average EOS wt with respect to wi, Awt = f{wi)Awi: 

f{wi) = VtQ[z + ?,wq (z - /n(l -h z)) (1 - Dq)] . (13) 

We consider wt because the measurements of are directly sensitive to 
Wt, so that models can only be distinguished if they have different wt- As 
can be seen from eq.(|T^), wt is much less sensitive to changes in wi when 
it is negative than when it is positive, mainly due to the value of Qq being 
larger for positive values of wi, We conclude that, in order to detect that 
Wq is time-dependent, it must be that the time-variation is large, roughly 
wi > 0.5, and it helps if wi is positive. This corresponds to the case where 
acceleration is becoming stronger as time evolves. 


6 Combining Supernovae with other Approaches 


Measurements to determine the EOS of the dark energy can be direct or 
indirect. Direct methods, such as SNIa observations, the Alcock-Paczynski 
(AP) test [^, and the cosmic microwave background (CMB) attempt to 
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measure the Hubble parameter H, its derivative H' and or some function 
of them. Indirect methods, such as structure formation aspects of the CMB 
and measurements of large scale structure (LSS) try to infer wq{z) from its 
effects on structure evolution. 

An example of a complementary observation is the the Alcock-Pazcynski 
(AP) test. The physical transverse size of an object is given by dr = dA^O = 
dA being the angular distance and A9 the observed angular size. The 
physical radial size is di? = / ^/^dr = For a population of 

spherical objects, the AP test is given by equating the transverse and radial 
sizes: AP(z) = ff = H(z)r(z) = H{z) f. 

The AP test on its own is not expected to improve the resolution of the 
dark EOS since it has a more complex dependence on wq than d^. What 
does seem promising, as pointed out by McDonald |^, is that the AP 
test can further constrain the range of Dm- 

Fig.0 shows the likelihood contours assuming optimistic anticipated er¬ 
rors over a continuous range between z = 0 and 2 ; = 2 of 1.4% for d^, and, 
for the AP test, 50 bins between redshift z = 1.5 and z = 3 measured with 
3% error per bin.|^ Both simulations represent highly optimistic assump¬ 
tions about future measurements. The results are interesting. The better 
constraint on from the AP test reduces the uncertainty in wq, but does 
not signihcantly change the uncertainty in the time-variation, Wi. This is not 
surprising since even a perfect determination of would leave a consider¬ 
able uncertainty in wi, as shown in Paper I. To be sure, the Alcock-Paczynski 
test is useful and worth pursuing, and a highly precise measurement com¬ 
bined with a highly precise measurement of SNe could determine the present 
value of w to within 15 or 20 percent. However, it does not help signifi¬ 
cantly with the particular problem of pinning down the time-variation of the 
equation-of-state. 

Measurements of the CMB anisotropy provide an additional probe of 
w{z). This probe also suffers from a degeneracy problem, even in the case 
where w is constant. The positions of the acoustic peaks in the temper¬ 
ature anisotropy power spectrum depend on the angular distance (d^) to 
the last scattering surface which, just like the luminosity distance for su¬ 
pernovae, depends on a multi-integral over w{z), dA = di/il + z)"^. In 
addition, the heights of the peaks depend on Dq and the Hubble parameter, 
Hq = h 100 km/parsec/sec. When all effects are considered, then, as shown 
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Figure 7: Two-sigma contours in the (tco,'W^i) = {wq{z = 0), dioq/dzo) 
plane for two idealized experiments. One measures thousands of supernovae 
between z = 0 and z = 2 (dashed contours). The supernovae are divided 
into 50 bins with a net error of 1.4% per bin. The second experiment is 
an optimistic estimate for the AP test (solid contours), assuming 50 bins 
of Lyman-alpha clouds uniformly distributed between z = 1.5 and z = 3 
with each bin measured with an accuracy of 3%. Both experiments assume 
a hducial model with = 0.3, Og = 0.7, wq = —0.7 = const., indicated by 
the X. In both experiments is marginalized over the range 0.2 to 0.4. The 
two-sigma joint likelihood for the two observations is shown in the shaded 
region. 
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by Huey, et al. |Q, the power spectrum is unchanged as certain combina¬ 
tions of Vtm, h, and w are varied. Consequently, none of these parameters 
can be determined well by the CMB data alone. Instead, measurements can 
only constrain these parameters to a thin two-dimensional surface in this 
three-dimensional parameter subspace. 

The reason why one might be optimistic about combining CMB anisotropy 
and SN measurements is that the degeneracy surface for the CMB anisotropy 
measurements is nearly orthogonal to the degeneracy surface for the SN mea¬ 
surements for the case of constant w. Figure 8 illustrates the small overlap 
between the SN and CMB degeneracy regions in the Vtm-w plane. Other au¬ 
thors have considered adding the CMB contribution!^, but they have not 
included the degeneracy aspect. As we shall show below, introducing time- 
varying w{z) introduces additional degeneracy that spoils the resolution even 
when the SN and CMB anisotropy measurements are combined. 

Rather than do another complete survey, which is a major technical chal¬ 
lenge on its own, we illustrate the degeneracy in parameter space with a 
simple example in which we consider the family of w{z) of the form: 


w{z) = Wq + WiZ foTZ<2 




(14) 


This form was chosen to allow signihcant time variation recently when Hg 
is large and, in particular, to be similar to the models considered in Paper I 
for < 2. For z < 2, the degeneracy problem with respect to SN data was 
already demonstrated and the wq-wi degeneracy region was characterized. 
However, we could not simply maintain the linear change in w{z) with respect 
to 2 : out to the last scattering surface at z = 1000 because the value of w 
would be ridiculously non-physical. Hence, we cutoff the ^(-dependence at a 
value of where fig is negligible and w{z) is physically plausible. We, then, 
maintain that condition back to the last scattering surface. For example, for 
w = —2/3 — l/6z and (Hg)o = .7, at z = 2 the dark energy contribution 
to the total energy density is less than 15%, which makes the details of the 
z-dependence cutoff unimportant. From z = 2 until last scattering surface, 
this model will have w = —1. 

The value of wq in our time-varying examples is hxed to be —2/3 except 
where otherwise stated. In each of these models, we also have h = .65, 
Qq = .7, Qm = -3, and Qf, = .04. Here Qb is the baryon density and Qm is 
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Figure 8: A simulation of the problem that arises if one assumes w{z) is 
constant in the htting procedure. For a given hducial model, the likelihood 
£t for the CMB anisotropy (dashed line) and SN luminosity distance-red shift 
(contour) observations are illustrated. The degeneracy curve for the CMB 
assumes cosmic variance limited sampling, and the SN contour assumes 1% 
error in luminosity distance. Each degeneracy region is long and thin, and 
the two are nearly orthogonal. Based on the small overlap, one is tempted 
to conclude that constancy of w is well established and its value is well 
determined. However, that conclusion is absolutely wrong. The hducial 
model in this example actually has a rapidly time-varying tn( 2 ;) = —2/3—l/6 z 
for z < 2 and w{z) = —1 for z > 2. The degeneracy regions were computed 
assuming tci = 0, but, if wi is hxed at a value somewhat less than zero, say, 
there are once again two narrow degeneracy regions which intersect over a 
small region, but the value of wq in the overlap region is significantly shifted. 
That is, the two experiments produce two degeneracy surfaces that intersect 
along a curve in the Wi direction along which a degeneracy remains. 
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CMB SNIa 



Figure 9: Illustration of the degeneracy problem for a model with constant 
w and two models with time-varying w as discussed in the text. The up¬ 
per left hand panel compares the CMB power spectra. The lower left shows 
the differences between the time-varying models and the constant w model 
and shows that they are less than or comparable to the full-sky cosmic vari¬ 
ance theoretical uncertainty, the envelope shown in the hgure (dotted lines). 
The upper right panel compares predictions for the luminosity distance-red 
shift relation. The lower right panel shows the differences with respect to 
constant model are less than the 1% resolution anticipated from supernovae 
measurements. 

the total matter density (baryonic plus non-baryonic). Note that luminosity 
distance-red shift measurements are not sensitive to Qb/^m, but the CMB 
measurements are. 

The time-varying models were treated as the hducial model, and then a 
numerical search was performed for a constant EOS model that is indistin¬ 
guishable from the hducial model based on the combined measurement of 
the CMB and of supernovae. Models were considered degenerate under the 
combined tests if: (1) the percent difference between the luminosity distance- 
redshift predictions for the two models is less than one percent out to z = 2 
(the same criterion as in Paper I); and, (2) the CMB predictions for the two 
models assuming a full-sky cosmic-variance limited measurement (no exper¬ 
imental error) cannot be distinguished to better than 3a. Both criteria are 
based on optimistic predictions of what will be realistically possible. 
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For the CMB, distinguishability between a model with a constant w and 
a hducial with a time-dependent w was determined by a log-likelihood analy¬ 
sis. The log-likelihood was calculated according to the log-likelihood formula 


obtained by Huey, et a/. [34 
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The coefficients C) and C) are the CMB multiple moments corresponding 
to the hducial and constant equation of state models, respectively. 

Fig. 8 illustrates the problem that arises if one assumes w{z) is constant in 
the htting procedure. We have already observed that this distorts results for 
the case of SN data alone. Here we show that the problem remains if CMB 
data is co-added. Assuming w is constant (wi = 0), both measurements pro¬ 
duce a thin degeneracy region in the plane. Based on the small overlap, 

one is tempted to conclude that constancy of w is well established and its 
value is well determined. However, this conclusion is absolutely wrong. In 
this example, the hducial model actually has a rapidly time-varying EOS 
w{z) = —2/3 — l/ 6 z for 2 ; < 2 (which produces a change in w of 50% over 
this range), and w{z) = —1 for z > 2. The degeneracy regions were com¬ 
puted assuming wi = 0. If wi were to be hxed at a diherent value, once again 
the two measurements will give two narrow degeneracy regions with a small 
overlap, but the value of wq in this overlap region is signihcantly shifted. For 
example, as shown in Fig 8 , hxing wi = 0 results in wq = —.74 with a few 
percent error. However, if Wi were to be hxed at its correct value Wi = —1/6, 
the result would have been wq = —2/3, again with a few percent error. But 
the central values diher by more than 10%. It is clear, then, that the two 
measurements produce two degeneracy surfaces which intersect along a de¬ 
generate curve which passes through a range of models with varying values 
for wq and Wi, that remain degenerate under the combined observations. For 
more complicated functional forms of wq{z) the degeneracy curve becomes 
a more complicated higher dimensional surface, and the range of degeneracy 
in parameter space (say, for wq) increases. 

The complications in the process of extracting the EOS from both mea¬ 
surements are further illustrated in Fig. 9. There we show two time-varying 
models with slopes [rcij > 0 . 1 , one of which is degenerate (by the log- 
likelihood test) with a constant w model with w = wq = —.72 and H;, = .04, 
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Figure 10: The same as Fig. 8 but with a hducial model with wi = 1/3. 
For cases like this with very rapid time-variation in ta, a symptom is that the 
CMB and SN degeneracy regions do not overlap. For wi large and positive, 
as in this example the SN contour (solid black) lies to the right of the very 
thin CMB degeneracy region (dashed curve). For wi large and negative, the 
SN contour lies to the left. 

VLm = .27, fig = .69, and h = .64. The other can be barely resolved making 
the most optimistic estimates about cosmic variance. A slight decrease in 
tci, or a slight decrease in experimental sensitivity would render the second 
model degenerate. The lower two plots magnify the differences between the 
predictions of the models. For the case of the CMB, we have also shown the 
envelope based on the constant w model corresponding to the full-sky cos¬ 
mic variance limit. For the SN, we have constrained the limits to he between 
± 1 %. 

If 1^11 is larger than 1/6 for our particular form of w{z)^ we hnd that there 
is no overlap between the degeneracy curve picked out by CMB measurements 
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and the degeneracy contour picked out by SN measurements (where both fits 
assuming w is constant). An example is shown in Fig. 10. In the case 
of negative (positive) Wi the CMB measurements that fit best suggest low 
(high) VLrn whereas the SN measurements suggest high (low) If this 

absence of overlap were to be found in the real data, an interpretation to 
pursue is that w is rapidly time-varying. Yet such an extreme scenario is not 
favoured by most theoretical models, most of which predict a moderately 
time-varying w. For the more likely case, in which the two measurements 
do overlap, combining them reduces degeneracy by only a modest amount, 
generally not even enough to decide whether the dark energy in the universe 
has a time varying equation of state or not. 

Co-adding the CMB to the SN data represents an improvement in the 
sense that Wi and Wq are more constrained than with SN data alone based 
on our earlier analysis or in Paper I. The improvement is by a factor of four 
or so assuming a linear form for w(z), which is signihcant. However, there 
remains a large uncertainty in the EOS. Furthermore, we would stress once 
again that the accuracy in determining w strongly depends on its assumed 
functional form. The range of degeneracy obtained for wi (a bit more than 
±0.1) in our example underestimates the degeneracy for general w(z). For 
example, for parabolic forms, the uncertainty in wi blows up to ±0.5. Given 
the extraordinarily precise data that has been brought to bear, the allowed 
variation in tco, and Qm is disappointing. 

7 Conclusions 

An important challenge for observational cosmology is to measure the equa¬ 
tion of state of the dark energy, wq(z). This can provide important infor¬ 
mation about the fundamental physics that is responsible for the accelerated 
expansion of the universe. Measurements of the distance-red shift relation 
using supernovae, perhaps combined with other direct methods such as the 
Alcock-Paczyinski test or the cosmic microwave background, would appear 
to be promising methods. Indeed, analyses based on the a priori assumption 
that wq{z) is constant suggest that wq can be resolved to 5% accuracy or 
better. 

In this paper and Paper I, though, we have uncovered a number of prob¬ 
lems and pitfalls that arise when trying to determine wq{z) without making 
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prior assumptions. Our lessons may be summarized as follows: 

• Because measures of luminosity or angular distance depend on inte¬ 
grals over wq{z), a first degeneracy problem arises in which neither the 
current value and nor its time-variation can be resolved to any useful 
accuracy. (Sec. II) 

• Since the effect of dark energy on the luminosity distance depends 
on the combination wqVLq rather than wq itself, a second degeneracy 
problem arises in which wq and Og are changed simultaneously so as 
to keep wqQq hxed. (Sec. Ill) 

• Although SN measurements may extend to z = 2, they are most sen¬ 
sitive to the behavior of wq{z) at a modest value of 2 :* ~ 0.1 — 0.4. 
(Sec. Ill) 

• Consequently, if there were only the hrst degeneracy problem, wq{z) 
could be well-resolved at z = z* even though it is not well-resolved for 
other values of ; 2 . Unfortunately, the resolution of wq{z*) is totally de¬ 
graded when one includes uncertainty in fig and the second degeneracy 
problem. (Sec. Ill, especially Fig. 2) 

• The common practice of htting data assuming that wq{z) is constant 
can lead to grossly distorted results. Similarly, the common practice of 
assuming wq > — 1 can lead to grossly distorted results. Fig. 5 shows a 
dramatic example in which these practices lead to the conclusion that 
Wq = — 1 and is well-resolved when, in reality, wq > — 1 and rapidly 
increasing. (Sec. IV) 

• Time-variation of wq is more easily detected if wq{z) is an increasing 
function of rather than decreasing. (Sec. V) 

• To resolve wq{z) with supernova data, an additional test is needed. 
Given optimistic estimates of experimental uncertainties, the Alcock- 
Paczynski test combined with the supernovae measurements can con¬ 
straint the current value of w to within 20 percent or so. However, Nei¬ 
ther the Alcock-Paczynski test nor microwave background anisotropy 
measurements provide the needed resolution to constrain the time- 
variation. (Sec. VI) 
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Our principal conclusion is that a new test is required to achieve the goal 
of measuring wq{z). In devising a new test, the two considerations must be 
precision and model dependence. Thus far, among the measurements that 
we have considered, the measurements which are precise give constraints on 
wq that are highly model dependent, leading to degeneracy problems. Tests 
which are not model dependent turn out to be difficult to measure precisely. 
So, there lies the challenge. 

In considering alternatives, it is critical to include practical estimates 
of their uncertainties. Furthermore, one must consider how the new tests, 
themselves, depend on wq{z). For example, claims have been made that Qm 
and Qq have been or will be measured very accurately by measurements of 
the cosmic microwave background [^. However, those estimates are based 
on assuming that wq = —1. Making no prior assumption about wq{z), a 
degeneracy problem once again arises]^ that spoils the resolution of Qq and 
wq{z), as discussed in Sec. VI. 

While trying to devise a new test to determine wq, it is worth mentioning 
that a precise measurement of H' will be extremely useful. The dependence 
of Wq on H, H' (prime denotes a derivative with respect to x = 1 + z) and VLm 

is given by wq = ^ 2 ^^ ^ good measurement of H' is clearly crucial 


0 


to the resolution of wq, but current tests do not probe H' directly. The 
next best option is to measure H{z), and then estimate H' by calculating its 
derivative. Obviously, this worsens the resolution for H' and increases the 
uncertainty in wq. 

Three additional approaches that we have not tried yet are measuring the 
time dependence of structure growth on gravitational lensing; and 

direct measurements of dz/dt (to be discussed elsewhere [p7|] ). 
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